On the validity of the adiabatic approximation in compact binary inspirals 
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Using a semi-analytical approach recently developed to model the tidal deformations of neutron 
stars in inspiralling compact binaries, we study the dynamical evolution of the tidal tensor, which 
we explicitly derive at second post-Newtonian order, and of the quadrupole tensor. Since we do 
not assume a priori that the quadrupole tensor is proportional to the tidal tensor, i.e. the so called 
"adiabatic approximation" , our approach enables us to establish to which extent such approximation 
is reliable. We find that the ratio between the quadrupole and tidal tensors (i.e., the Love number) 
increases as the inspiral progresses, but this phenomenon only marginally affects the emitted gravita- 
tional waveform. We estimate the frequency range in which the tidal component of the gravitational 
signal is well described using the stationary phase approximation at next-to-leading post-Newtonian 
order, comparing different contributions to the tidal phase. We also derive a semi-analytical expres- 
sion for the Love number, which reproduces within a few percentage points the results obtained so 
far by numerical integrations of the relativistic equations of stellar perturbations. 

PACS numbers: 04.30.-w, 04.25.Nx, 04.25.dk 



I. INTRODUCTION 

Coalescing binary systems of neutron stars (NS) 
and/or black holes (BH) are among the most interest- 
ing sources of gravitational waves (GWs) to be detected 
by advanced Virgo and LIGO One of the key fea- 
tures of the coalescence is the NS tidal deformation, 
vifhich provides precious information on the NS equa- 
tion of state (EOS) 0-0] ■ For this reason, theoretical 
and numerical studies have been recently performed to 
model the effect of tidal deformations on the emitted 
GW signal, and to extract its contribution from a de- 
tected signal [1,11, [B^ni- These studies are based on fully 
general-relativistic numerical simulations and on (semi- 
)analytical approaches. 

In current analytical approaches, the tidal deformation 
properties of NSs are encoded in a set of numbers, the 
Love numbers [Hi Ell i which relate the mass mul- 
tipole moments of the star to the (external) tidal field 
multipole moments. In the so-called "adiabatic approxi- 
mation" , at the lowest multipole order, the evolution of 
a star in response to an external quadrupolar tidal field 
Cij , is governed by the equation 



Qij — , 



(1) 



where Qij is the star quadrupole moment, i?NS is its ra- 
dius at isolation, and k2 is the (dimensionless) I = 2 apsi- 
dal constant, also dubbed second tidal Love number. For 
non-rotating relativistic stars, k2 is usually computed by 
solving the linear I = 2 static, even-parity perturbations 



of Tolman-Oppenheimer-Volkoff solutions |14| . requiring 



regularity and continuity for the metric perturbation and 
its first derivative 

IlMlil- Eq. dH) and its higher mul- 
tipole order versions have been employed to determine 
the effect of tidal deformations on the orbital motion, 
and on the GW signal emitted by NS-NS and BH-NS bi- 
nary systems [1, (a^H, |1] . In these studies, the relation in 



Eq. ([T]) — or an equivalent assumption [sj — is assumed 
as a starting point. 

The domain of validity of ([T|) is briefly discussed in 
, using previous results of Lai. In [11] a NS tidally in- 
teracting with a companion, is described as a forced os- 
cillator in a Newtonian framework; the energy absorbed 
by the oscillator is the sum of an "instantaneous" term, 
proportional to the forcing (tidal) field, and a term as- 
sociated with the stellar oscillations. Since the former 
is much larger than the latter, it is argued that the adi- 
abatic approximation ffl holds [1]. A forced oscillator 
model is also used in [3|, providing further evidence of 
the accuracy of the adiabatic approximation. It should 
be stressed, however, that only a consistent dynamical 
study of the stellar deformation during the inspiral may 
assess the validity of the adiabatic approximation. This 
is the approach we follow in this paper. 

In a recent work [ll], some of us developed a semi- 
analytical description of tidal deformations of NSs in in- 
spiralling BH-NS and NS-NS binaries. This model com- 
bines the post-Newtonian (PN) approach (see, for in- 
stance, |17Hl9f). which accurately describes the orbital 
motion before the onset of mass-shedding, and the ajfine 
model, which allows for a description of stellar deforma- 
tions due to an external quadrupolar tidal field generated 
by the companion j2(]| - |25| . We computed the tidal tensor 
associated with the PN metric of a two-body system, de- 
fined in terms of the PN Riemann tensor and of the NS 
local tetrad 



T-) a l3 ^ 5 



(2) 



up to 1.5PN order. This tensor appears as a source in 
the dynamical equations describing the stellar deforma- 
tions. Note that the use of a two-body PN metric explic- 
itly yields self-interaction terms in the tidal tensor. We 
validated our post-Newtonian-affine (PNA) approach by 
comparing the results obtained for BH-NS binaries to 
the outcome of fully general-relativistic numerical simu- 
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lations. In addition, we computed the Love number k2 
using Eq. ([1]) at large separations, showing that our re- 
sults were in good agreement with the analytical results 
of [11. 

In this paper, we extend the PNA approach further: 

• We compute the tidal tensor ^ up to 2PN order. 
This new result may allow one to construct more 
accurate models of gravitational waveforms emitted 
by inspir ailing compact binaries, which could be 
employed to extract information on the Love num- 
ber and on the underlying equation of state, dur- 
ing post-processing. Indeed, the 2PN terms may 
be used to correct the binary binding energy and 
the GW-flux and to improve the description of the 
phase evolution. 

• We use the 2PN tidal tensor to assess the range of 
validity of the adiabatic approximation. We do so 
by solving the dynamical equations for the orbital 
motion and the stellar deformations, by computing 
the quadrupole and tidal tensors, and by determin- 
ing /c2 from 



3Q. 



(3) 



at different values of the orbital separation. We 
find that k2 increases in the late inspiral phase; it 
should thus be referred to as a Love function k2 (r) . 
of which k2 is the asymptotic limit. As discussed in 
Section IIIIl we prefer to express /c2 as a function of 
r, rather than as a function of the gauge invariant 
frequency, because in this way it is easier to find 
an accurate fit for k2. We then use the station- 
ary phase approximation [2^ |23 to compute the 
gravitational waveforms with tidal phase effects in- 
cluded up to IPN relative order Q and model the 
gravitational wave phase accordingly. 

• We compute the fitting factors between point- 
particle templates and a model of the gravitational 
signal that includes tidal effects at the best of our 
knowledge up to IPN order, i.e., including the 
Love function k2(r). We find that point particle 
templates satisfy the accuracy standards defined 
in [2^, unless one considers NS-NS binaries with 
a very stiff equation of state, in which case the 
aforementioned standards are marginally violated. 
Moreover, we show that tidal effects can affect the 
measurement of the total mass and the symmetric 
mass ratio by at most 3% and 2%, respectively. 

• We derive a simple, semi-analytical expression for 
the tidal Love number k2'. 



k2 
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of the pressure over the stellar volume, both cal- 
culated for the star at isolation. We check the 
accuracy of this formula by comparing its results 
to those obtained by perturbative approaches [l^ . 
finding that they agree within a few percentage 
points. 

The plan of the paper is the following. In Sec |ll] we 
briefly describe our model. In Sec. IIIIl we study the tidal 
Love number k2 in the post-Newtonian affine approach. 
In Sec. llVl we compute the gravitational waveform in the 
stationary phase approximation, comparing the different 
tidal contributions to the Fourier phase. In Sec. |V] we 
compare the waveforms obtained by using different ap- 
proximations for the tidal contributions to the Fourier 
phase. In Sec. IVII we draw the conclusions. 



II. THE MODEL 

In the following we briefly sketch the post-Newtonian- 
affine model. For further details, see [1^1 ■ The masses 
of the two compact objects, inspiralling on quasi-circular 
orbits, are mi and m2, m = mi + TO2 is the total mass, 
and u = mim2/m? is the symmetric mass ratio. We 
write the equations for the secondary object m2, of radius 
i?NS- In a NS-NS system, analogous equations hold for 
the primary object mi. 



A. Tidal deformations in the affine model 

To describe the stellar deformation, we use the aflane 
model approach H^-Hl (improved in jll,[2l[2l|), which 
is based on the assumption that a NS in a binary system 
preserves an ellipsoidal shape when it is deformed by the 
tidal field of the companion. The deformation equations 
are written in the principal frame, i.e., the one comoving 
with the star and with axes coincident with the principal 
axes of the ellipsoid. Surfaces of constant density inside 
the star form self-similar ellipsoids, and the velocity of a 
fiuid element is a linear function of the coordinates cc* of 
the principal frame. Under these assumptions, the infi- 
nite degrees of freedom of the stellar fluid reduce to five 
dynamical variables and we are left with a set of ordinary 
differential equations describing the evolution of the star. 
The five variables describing the stellar deformation are 
the three principal axes of the ellipsoitfl fli (* — 1,2,3) 
and two angles ip, A defined as 



(IT CLT 



(4) 



where M. and 11 are, respectively, the scalar 
quadrupole moment of the star and the integral 



^ In our conventions, '1' denotes the direction along the axis di- 
rected towards the companion, '2' is the direction along the other 
axis that lies in the orbital plane, and '3' indicates the direction 
of the axis orthogonal to the orbital plane. 
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where r is the NS proper time, is the elhpsoid angular 
velocity measured in the parallel transported frame asso- 
ciated with the star center of mass (the "figure" velocity 
in [l^l), and A (defined in [22[ in terms of the vorticity 
along the z— axis in the corotating frame) describes the 
internal fluid motion in the principal frame. 

The NS internal dynamics is described in terms of the 
Lagrangian 



Ti-U -V 



(5) 



where Tj is the star kinetic energy, U is the internal en- 
ergy of the stellar fluid, and V is the star self-gravity. 
The relevant integrals are the scalar quadrupolc mo- 
ment Ai = 1/3 / (iMB^j(xi)^, the pressure integral 
n = JdMsp/p, and the self-gravity potential V = 
— / dMs J2i^i^i^' here dMs is the baryon mass ele- 
ment, $ is the gravitational potential, p is the pressure, 
and p is the baryon mass density. 

As shown in [l^, [2l| , under the affine hypothesis the 
integrals , 11, V can be expressed in terms of integrals 
over the spherical configuration of the star (a^ = i?Ns)j 
and of functions of the dynamical variables. In the fol- 
lowing, all carets ( " ) denote quantities computed on the 
spherical star. In the spherical configuration, the pres- 
sure integral and the self-gravity are related by the virial 
theorem 



V -3n. 



(6) 



The original afhne app roach, introduced in a Newto- 
nian framework [20|-|22|, was improved and extended in 
(igI [23 . [25j |. The spherical configuration of the star is 
determined by solving the rclativistic equations of stel- 
lar structure, which yield the profile of p(rs), being rg 
the radial coordinate in a Schwarzschild frame associ- 
ated to the NS. The gravitational potential $ appearing 
in the definition of V is replaced by an effective rcla- 
tivistic potential such that the virial theorem ^ is sat- 
isfied. The equations of motion for the internal vari- 
ables qt = {'0, A, ai, 02, as} and their conjugate momenta 
Pz = {Pi^,P\,Pai,Pa2,Pa3} are: 



da I 
~dt 


_ RmPai 
7(i) M 


(7) 


da2 


RnS Pa2 
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dpai 

dt 



dPa^ 

dt 
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In the above equations, the c^'s are the PN tidal tensor 
components defined later in Sec. Ill CI 



du 



lo iaf + u)y/{al + u){a2 + u)(a| + u) 
are elliptic integrals. 



(17) 
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n = ^4. 
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(20) 



where xa = ruA/m {A = 1,2), \/(r,) ^ G^^ ^dr^, 

msirs) is the gravitational mass enclosed in a sphere of 
radius , r and v are the orbital distance and the relative 
velocity of the two bodies in the PN frame, and 



C 



J2 



M 

p2 

M 



[[af + al) A - 2aia2r2] 
[{al + al)fl- 2aia2A] 



(21) 



are the circulation and the NS spin angular momentum. 
In this article we consider the case of an asymptotically 
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non-rotating star. As discussed in [23 |. the intrinsic spin 
of the star may be defined only far away from the com- 
panion, where the star is axisymmetric and the elhpsoid 
rotation n vanishes; in this region, it is possible to iden- 
tify the angular velocity of the star with —A. Therefore, 
an asymptotically non-rotating star has A = 0, = 0, 
and thus C = 0, at r — >■ oo. Since C is a constant of 
motion, it remains zero during the inspiral. 51 and A, in- 
stead, become non- vanishing as the NS couples with the 
tidal field, i.e., the ellipsoidal star acquires a (very small) 
angular velocity associated with general relativistic ef- 
fects, such as geodetic precession and frame dragging. 



B. The orbital motion 

The orbital motion is described in the post-Newtonian 
framework starting from the the 3PN metric of the two- 
body system written in harmonic coordinates {x" = 
ct,x,y,z} [la, We assume an adiabatic inspiral of 
quasi-circular orbits, i.e. such that the energy carried 
out by GWs is balanced by the change of the total bind- 
ing energy of the system [l^. We use the Taylor T4 ap- 
proximant to describe the orbital phase evolution of the 
two-body system, including the effects of the tidal inter- 
action on the orbital motion, derived at the IPN order 
(beyond the leading term) in The phase 0(t) and the 
orbital frequency uj = d(j){t)/dt are found by numerically 
integrating the following ordinary differential equations: 



dx _ 



pp 



f -^tid 



dt Gm 



where 



/ Gmuj\ 



2/3 



(22) 
(23) 

(24) 



termine the radial coordinate r{t), we employ the PN 
expression for 7 = (not to be confused with the time 
dilation factor 7(i) defined in Eq. ((T8|) ). which is also 
known up to 3PN order [s^l- We refer to [l^ for the 
explicit form of the equations, or to the original papers 



C. The post-Newtonian tidal tensor 



The quadrupolar tidal tensor is 



(0)e( 



^(o)^(i) 



(25) 



where Rap-yS is the Riemann tensor of the 3PN metric de- 
scribing the orbital motion. We introduce an orthonor- 
mal tetrad field e(i) {i = 0, . . . , 3) associated with the 
frame fixed to the star center of mass O* , parallel trans- 
ported along its motion, and such that e^^-j ~ w^, i.e., 
the 4— velocity of O* . In [l^ we computed Cij up to 
0(1/0'^), whereas previous computations were performed 
up to 0(1/0^) [13, [1H|- In this work we further expand 
the PN tidal tensor in order to include also 1/c^ terms. 
To this aim we first need to calculate the 2PN component 
of the spatial tetrad vectors e-;!^.-, introduced in [1^ . From 

Sij we find that 
the tetrad field expanded up to 1/ order for the bodies 
A = 1,2 is: 



(fc) 

the orthogonality condition g^j/e^^^e^^-j 



Tpp is the point-particle term [32| , and J-tid incorporates 

the finite-size effects on the orbital motion [^. To de- with 
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(26) 



-A{V,)A + v\['mA + ^ 



{V)a 



^A^A 

2c2 



+ 0(c-5) 

mlSl + iv\v\ {A{V)a + v\) - lQ{W,k)A] + 0{c-^) 
I 



where {V)a, {Vi)A, {Wij)A are PN potentials defined in the body A, and regularized as described in [T^; finally, 
terms of the source densities, evaluated at the location of 



5 



is the angle describing geodesic precession and frame with to arbitrary initial time. 



dragging 



(Q 



Following the same procedure described in [16| we 
project the Riemann tensor on the tetrad field, to de- 
rive the tidal tensor which we express in the NS principal 
(27) frame as c = TCT^, where 



here Qxy is the non vanishing component of the antisym- 
metric matrix Q defined as 361 : 



cos ip sin 'ip 
T — \ — sin Ip cos Ip 
1 



(29) 



Q{t, to) = / [v X (Vy - a) - V X (Fv - 2V)] dt , 

J to 



and ip has been defined in Section Hi Al We find that the 
tidal tensor components for the body 2, appearing on the 
(28) right-hand-side of the dynamical equations pUjl - p^ . are: 



^ {1 + 3cos[2^,]} + ^ {^(5 + xO (1 + 3cos[2^.]) - 6v^ (1 + cos[2^,])} + {^X + 
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. 2 / 223 3 , 
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-(xi-5x2)cos[2V.,]} + ^. 
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+ ^ [2X2 - 2x2 + 1] J> sin[2^,] -f 0{c-') . (33) 



In the above expressions, 6x = Xi ^ X2 and the lag angle 
ipi = Ip — 4> + ^ describes the misalignment between the 
ai axis and the line between the two objects. The super- 
script dot identifies derivatives of the orbital variables (p 
with respect to the coordinate time t. It should be noted 
that in the center-of-mass frame of the binary system, the 
tidal tensor is given by the same expressions ([50]) - ([55)1 . 
setting ipi = —(p. We also stress that, by construction, 
the tidal tensor is traceless. 



Finally, we remark that (for a star with zero circula- 
tion, i.e., C = 0) the lag angle is very small, since it is 
only due to general relativistic effects such as geodetic 
precession and frame dragging; therefore, '0; <C 1. Then, 
since at leading order Cxy oc sinp?/)/], it follows that 

We have checked numerically that Cxy is always several 
orders of magnitude smaller than the others components 
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of the tidal tensor. We remark that, as noted in [37| . 
the lag angle becomes non-negligible when viscosity is 
included in the model. 



III. TIDAL LOVE NUMBERS IN THE 
POST-NEWTONIAN-AFFINE APPROACH 

NS tidal deformations can be described in terms of a 
set of parameters, the Love numbers, encoding the defor- 
mation properties of the star. The idea is that in pres- 
ence of a weak, external tidal field, a spherical star is de- 
formed, and its mass multipole moments are proportional 
to the multipole moments of the perturbing tidal field. In 
this paper, we focus on the Love number k2, which is as- 
sociated with the lowest order "electric" moment, i.e., 
the ^ = 2, or quadrupole, moment. This is the most rel- 
evant for the phenomenology of stellar deformations, to 
be considered when modelling gravitational waveforms. 
The mass quadrupole moment (traceless) tensor is pro- 
portional to the tidal field 0, : 

Q.J = -^fei?^sQ, . (35) 

In the NS principal frame, this equation reads 

2 

qtj = -gfc2i?NsCjj , (36) 

where qij is the quadrupole moment tensor projected 
onto this frame 

<?y = J dAlB (^x'x^ - ^S'^r^^ d^x 

= ^i<A-^'W' (37) 

(no sum on i), — {a\ + a\ + a3)/3, and M is the 
scalar quadrupole moment for the spherical, isolated star 
configuration given in Eq. (|19|1 . 

Eq. p6p is based on the adiabatic approximation, 
which assumes that the timescale of the orbital evolu- 
tion (and thus of tidal tensor variations) is much larger 
then the timescale needed for the star to readjust into a 
stationary configuration. In the following we will drop 
this assumption and assess its validity. 



high compactness C (not to be confused with the circula- 
tion defined in Eg. 121(1 and small deformability A, whereas 
the PS EOS describes stiff NS matter and gives stellar 
models with small compactness and large deformability 
(see for instance [^). We consider two values of the sec- 
ondary NS mass, m2 = 1.2 A/q, 1.4 M0, and three values 
of the mass ratio, q = mi/7712 = 1,3,5. In the q = 1 
cases, both bodies are NSs, with the same EOS, and we 
take into account the deformation of both of them. When 
(7 = 3, 5 we assume that primary body is a BH. 

For each of these binary models, we compute k2 from 
Eq. dSnH, using Eq. ^ for the left-hand side. We find 
that this quantity is a function of the orbital distance, 
k2 (r) , and that it increases during the stellar inspiral. We 
express k2 as a function of r, and not of the PN variable x, 
because it is easier to find an accurate fit for ^2 in terms 
of the r-dependence rather than of the x-dependence. 
Further, we think that the variable r provides a more di- 
rect physical insight on the process. Hereafter, we denote 
the asymptotic limit of the function k2{r), i.e., the Love 
number, as 

k2 = hm k2 (r) . 

J — >oo 

We determine k2 by fitting our data with the following 
expression: 

k2{r)^k2{l + a-+p—+^—]. (38) 

In Fig. [1] we plot the function ^2 {r) for all the binary 
configurations we considered. The figure shows that in 
the last stages of the inspiral, before the mass-shedding 
sets in, k2{r) grows by a factor of ~ 10% — 30%. This 
effect is strongly dependent on the EOS choice: it is 
larger/smaller for a stiffer/softer EOS as PS/ APR, or, 
equivalently, for a less/more compact NS. The depen- 
dence of this effect on the mass ratio q is weaker; the 
Love number increases more for larger values of the mass 
ratio, unless the star reaches the innermost circular or- 
bit (ICO) before the mass-shedding sets in. It seems 
that this effect is nearly insensitive to the NS mass. We 
show the curves in Fig. [T] as functions of r, and not of 
r/m. because we find that this choice makes the common 
asymptotic limit of the curves more evident. 



A. Evolution of the tidal Love number during the 
inspiral 

We solve Eqns. dTD-dH]), and (gSl), for a represen- 
tative set of binary system configurations. We consider 
two equations of state, dubbed APR (which had been 
derived by Akmal, Pandharipande and Ravenhall [37]) 
and PS (which had been derived by Pandharipande and 
Smith [38]). These EOSs are expected to cover the range 
of possible EOS stiffness and NS deformability: the APR 
EOS describes soft NS matter and yields NS models with 



B. Equations of tidal deformation in the 
perturbative regime 

In the limit of small deformations 

Sal = al~ i?2 s < , (39) 

it is possible to derive a semi-analytical, closed form for 
the Love number ^2 , expressed in terms of the NS scalar 
quadrupole moment and of the integral of the pressure 
over the stellar volume. 



7 



APR, m2=1.2Mg 
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FIG. 1. (Color online) Love function k2{r) for different choices of the NS EOS, NS mass m2, and binary mass ratio q, as 
indicated in the panels. For the APR EOS, q = 5, and both values of the NS mass, the NS reaches the innermost circular orbit 
before the mass-shedding sets in. 



By introducing the notation Sa^ = ^ ^'^i 



i?^g, wc rewrite Eq. ([571) 



M 



{5a1 - 50^)5'- 



(40) 



(no sum on i). As a first step, we show that the spherical, 
non rotating solution = i?NS- = A = 0, is a solution 
of Eqns. ([71)- P^ . with Cij = 0. Indeed, in this limit the 
dynamical equations reduce to 



M 

lit) 



1 V 



27W 



—_R%aA 



n 

7^ 



where we introduced A = Ai{ai 
Eq. ((ni) thus reads 



= , (41) 
^Ns), which from 



A^ 



du 



2 1 



(i?2 



NS 



(42) 



Eq. dH]) is therefore satisfied when Eq. (O, V = -3n, 
i.e., the virial theorem, holds. 



We now consider the first order perturbative expan- 
sion of Eqns. ([7)l- (|16p around the spherical, non- rotating 
configuration. This expansion is an accurate description 
of the dynamical system when deformations are small, 
i.e., when the tidal field is weak and the rotation rate is 
small. 

To simplify the discussion, we restrict the discussion 
to the case of an asymptotically non-rotating star; in 



and 



this case Eq. ([M)l holds, i.e., Cxy <C „„, 
Eqns. (HHl-llTll) and (HI]) guarantee that Vl, A, A, and V 
may be neglected. The remaining non-trivial equations, 
([Tni)-(|I21), reduce to 



2 M 



i?lg II 

M af 



= 



(43) 



{i = 1, 2, 3, no sum on i). We remark that in this proof 
the explicit expression of the tidal field c^, which drives 
the stellar deformation, can be taken as generic (we only 
require that c^y < Cxx,Cyy,Czz)- 
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Expanding Eq. ([^S]) wc find 







2Sa^\ 




2M 




5 ^NS / 


M \ 



7?2 



(44) 

where we denoted the expansion of 11 as 11 + SH and 
considered the Taylor expansion of the integrals Ai P7|) 

A = ,45) 



5 i?^g 



Imposing the virial theorem V = —311, one now has 
n /SU 2 Sa^ 



n 

The traceless part of Eq. (|46)) yields 
2 ft (5a? - (5a2 



0. 



(46) 



i?2 



and plugging this into Eq. (|40[) gives 



(47) 



(48) 



A similar relation was found in [29[ . 

Finally, Eqns. (|36)l and (j48p yield the semi-analytical 
expression for the Love number we were looking for: 
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(49) 



where fl and M are given in Eqns. ([^0]) and (|19p . respec- 
tively. 

In Table |T] we show the quantities characterizing the 
different binary models considered in this paper, and 
the corresponding values of the Love number ^2, com- 
puted in three different ways: (i) from the dynamical 
evolution of our equations ©-(Ull), (ii) by 

solving the equations of relativistic stellar perturbations 
for an isolated NS, derived in [T^] (^2^); (iii) from the 
semi-analytical forrnula ((49|) (fcf")- We notice that, as 
expected (see also [IGj , where polytropic EOSs were con- 
sidered), the less relativistic the NS is, i.e. the lower 
its compactness is, the lower the relative error between 
the relativistic value and our perturbative result fc|" 
is. Note that when the Love number is extracted from 
the dynamical evolution of the binary system, it has a 
(weak) dependence on the mass ratio q; the quantities 
k2 , fcf"; instead, do not depend on q, because they are 
evaluated in terms of the intrinsic properties of the star. 
We find that the values of k2 computed with these three 
approaches have very small discrepancies, of at most few 
percentage points. 

Table H] also provides the value of the orbital frequency 
fcut where our simulations stop. This corresponds to 
the onset of mass-shedding, or to the ICO, if the latter 
is encountered before mass-shedding (see discussion in 

M). 



EOS m2(M0) 


q 




Afc2/fe2 






APR 


1.2 


1 


0.0884 


0.011 


0.031 


532.43 




1.2 


3 


0.089 


0.019 




500.32 




1.2 


5 


0.090 


0.033 




404.88''' 




1.4 


1 


0.079 


0.047 


0.058 


575.88 




1.4 


3 


0.081 


0.071 




541.07 




1.4 


5 


0.080 


0.063 




347.O7I' 


PS 


1.2 


1 


0.117 


0.024 


0.006 


354.81 




1.2 


3 


0.114 


0.048 




339 3Q 




1.2 


5 


0.116 


0.034 




325.55 




1.4 


1 


0.110 


0.017 


0.017 


378.04 




1.4 


3 


0.108 


0.032 




354.11 




1.4 


5 


0.112 


0.005 




345.20 



TABLE I. Values of the Love number for different EOSs, NS 
masses m2, and binary mass ratios q (columns 1,2,3, respec- 
tively) . The value of ^2 given in column 4 is computed solving 
Eqns. ©-(Uni), (1221), and ^E^; in columns 5 and 6 we show 
the relative errors {k2 — k^)/k2 and (fcf" — ^I')/^, where ^2^ is 
the relativistic Love number computed solving the relativistic 
equations of stellar perturbations 12], and is evaluated 
using the semi-analytic formula H49[) (see text). The values of 
the orbital frequency f°^t at which our simulations end are 
also provided in the last column. The ^ symbol indicates that 
the ICO is reached before the onset of the mass-shedding. 



IV. GRAVITATIONAL WAVEFORM IN THE 
POST-NEWTONIAN-AFFINE APPROACH 

We now compute the gravitational waveform including 
tidal effects by means of the Love function k2{r). In the 
following derivation, we use geometric units G ~ c = 1. 
Our starting point is the state-of-the-art inclusion of tidal 
effects in the GW signal by means of the Love number 
^2. This is based on the binding energy £{x) and the 
gravitational flux £gvf (a;) 0, HH: 



£{x) = - 



fe/2 T Al X 
— A2 

A2 m 




2 — — X 



+ 1^2} ,(50) 



5 X2m5 



6(3-2x2) 



1 



(704 + 1803x2 



28 

4501x2 + 2170x2) x]+1^2} , (51) 



where Xa (^ = 1,2) is the tidal deformability associated 
with the A-th body and related to the Love number ^2 A 
by 



Xa 



:h,AR'' 



NS 



(52) 
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The point-particle coefficients , are given in tlie Ap- 
pendix of |32| . 

Using tlie stationary phase approximation (SPA) (26l . 
[27| and the TaylorF2 framework to construct the 
signal in the frequency domain, the {I, mi) mode of the 
GW signal reads 



271 



mi(j){tf) 



(53) 



where i/ is defined as the time when the instantaneous 
frequency matches the Fourier variable, i.e. miuj(tf) = 
27r/, is the orbital phase, the dots indicate a second 
derivative with respect to time, A''"^' is the Fourier am- 
plitude, and -0'™' = 27r/i/ — 'ini<f>{tf) — j is the Fourier 
phase. To compute the (/) = ui entering the GW signal, 
we follow the strategy adopted in [sl] and express uj as 
u! ~ 2m V^-'^ J with X derived using the TaylorT4 prescrip- 
tion given by Eq. ([^D). Hereafter, we shall focus on the 
I ~ nil = 2 mode of the gravitational waveform and drop 
the superscript "22" for sake of simplicity. The coeffi- 
cients of the 3PN order expansion of the amplitude A{x) 
are collected in (s^]. The Fourier phase has the form 



"0 = V-'PP + 0T , 



(54) 



where the coefficients of the 3.5PN order expansion of 
the point-particle contribution i/)pp(x) are given in [32l | 
and the tidal term, -^t, is calculated up to IPN (relative) 
order in 0, H, |^, by assuming a constant Love number. 
It reads 



(55) 



with 



V'i™(/) = - 



3x5/2 
128i^to5 

3:^5/2 
128i/m5 



A, 



A, 



^(1 

X2 



llXi 



1 2, 



28X2 



(3179-919x2 



-228x^ + 260x^) x\+1^2. 



(56) 



(57) 



We may now use the Love function we introduced in 
Sec. mil as an "upgrade" of the concept of Love number 
and determine new tidal contributions to the GW signal 
emitted by an inspiralling binary. Its r-dependency in 
Eq. ([55)1 may be cast into a dependency on the PN di- 
mensionless variable x by using the PN expansion of m/r 
in terms of x j33l | . This way one has 



k2{x) = h [1 + ax] + Oix^) , 
and consequently 

X{x) = X[l + ax]+0{x^), 



(58) 



(59) 



where the coefficient a is determined by means of our 
dynamical simulations. The expanded AA(a;)'s may now 



replace the A^^'s in the PN expanded energy S{x) and 
GW flux Cgw{x), Eqs. dSOD and ^ respectively. No- 
tice that since the tidal correction to the binding energy 
and gravitational flux are known up to the IPN order. 
Eg. ([55]) was truncated at 0{m/r) = 0{x). We remark 
that the GW flux depends on the time derivatives of the 
quadrupole tensor of the system and, therefore, the 
replacement A^ —J- Xa{x) in Eqns. ([50]) and (|5ip neglects 
terms arising from time derivatives of Xa{x). Such terms, 
however, being proportional to the velocities of the com- 
pact objects, are of a PN order higher than the order of 
the expansion in Ea. (j57p . and may be safely neglected. 
The replacement Xa — >■ A^ (x) yields the following correc- 
tion to the Fourier phase of the gravitational waveform 
given in Eq. ([54]) : 



3x5/2 

128i/m5 



A2 5^(2 + 27— ]a2X +lo2, 

7 V X2/ 

(60) 

where aA are the coefHcients appearing in Eq. ([55]) asso- 
ciated to the A — th body and Xa ai'c determined from 
Eq.(l52]). Eq. dMl) then becomes 



27— 

X2 



^ = 0PP + ^ ^1,PN , 



(61) 



We remark that the term il))^^ and our correction (5-0T 
have different origin. The former, obtained indepen- 
dently in and in [U is the next-to- leading or- 
der tidal correction to the Fourier phase and was de- 
rived under the adiabatic approximation Qij = XCij (or, 
equivalently, assuming that the contribution of "electric" 
type quadrupole deformations to the action is AS* ^ 
XCijCij); the latter is the result of a dynamical evolu- 
tion in which we monitor the variation of the ratio be- 
tween the quadrupole and tidal tensor during the inspiral 
(which we find to be the same, within < 1%. for all com- 
ponents ij); this variation is encoded in the parameter a 
appearing in Eq. (j59p . Thus, these terms (both increas- 
ing as r decreases) simply add linearly. 

In Fig.[5]we compare the different tidal contributions to 
the GW Fourier phase: the leading order term (solid 
black line); the next-to- leading order 0t = V't + iPt"^ 
(dotted blue line); the total tidal contribution including 
our correction i/'T+i^V'T, (dashed-dotted red line) . As dis- 
cussed above, for consistency reasons our approach only 
allows us to include 1/r (^ x) terms in the gravitational 
waveform, therefore we truncate our expansion at order 
0{x). 

For each model the phases are shown in the range 
[10 Hz, /cut], where /cut = 2/°^^, (^^^ Table©. Fig. [H 
shows that: 

• The dcphasing due to the NS tidal deformation is 
significantly larger than one radiant only for NS- 
NS binaries. For BH-NS binaries, V't ^1 rad and 
therefore, as previously established in [6[, it is un- 
likely that second generation detectors such as Ad- 
vanced LIGO /Virgo will detect tidal deformation 
effects in GWs emitted by mixed binaries. Only 
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APR, m2=1.4Mgun, q=1 

4.5 r 1 1 1 r 




200 400 600 800 1000 

fGw[Hz] 



PS, m2=1.4M3,„, q=1 

12 r 1 1 1 1 ^ 




100 200 300 400 500 600 700 

fowlHz] 



APR, m2=1.4M3yn, q=3 



PS, m2=1.4M3,„ q=3 




200 400 600 

fGw[Hz] 



300 400 500 600 

fGw[Hz] 



APR, mo=1.4M, 



PS, m2=1.4M3,„ q=5 




300 400 

fGw[Hz] 



300 400 

fGW [Hz] 



FIG. 2. (Color online) Tidal contribution to the Fourier phase of the gravitational wave signal, for m2 — 1.4 Mq, the APR/PS 
EOS (left/right panels), and, from top to bottom, q — I, 3, and 5. Different curves show: the leading order term i/j^ (solid 
black line); the next-to- leading order contribution tpT = V't +V't^^ (dotted blue line); the total tidal contribution including our 
correction ipT + S^t, (dashed-dotted, red line). For each model the phases are shown in the GW frequency range [10 Hz, /cut], 
where /cut = 2/cut , /cut" being the orbital frequencies given in Table U 
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third generation detectors like the Einstein Tele- 
scope (ET) [4^, which could detect signals with 
very large signal-to-noise ratio, would be able to 
reveal the tidal contribution in BH-NS binary coa- 
lescences. 

• Looking at systems with q = 1, i.e., NS-NS bina- 
ries, the tidal signal is larger for stifFer EOS, or 
equivalently less compact NSs, and for lower values 
of the NS mass, confirming previous results [3ll4^. 

• The correction derived in this paper, SipT, affects 
the GW phase only marginally for the binaries we 
consider. Indeed, the two curves describing tpT and 
ipT + SipT nearly coincide. 

Similar conclusions hold when we choose m2 = 1.2Mq. 

We remark that in our approach we determine the 
gravitational waveform up to the frequency /cut, corre- 
sponding to the onset of the mass-shedding, which (in our 
PNA approach [l^) occurs when the NS fills its Roche 
lobe. As shown in Table HI we find values of /cut between 
^ 700 and ^ 1150 Hz. These values are larger than 
the cutoff frequency assumed in Q (/cut = 450Hz), but 
more conservative than the values employed in j4^ , which 
correspond to the configuration in which the surfaces of 
the two bodies touch. As discussed in Q, the main rea- 
son behind the /cut = 450Hz choice is that the next-to- 
leading order (IPN) corrections to the tidal phase, which 
we include, were neglected. On the other hand, the large 
values for /cut used in (see also [ll|) correspond to 
the configuration in which the surfaces of the two bodies 
touch. We think that ending the integration when the 
NS fills its Roche lobe is a safer choice. 



V. WAVEFORM COMPARISONS 

We now examine the dephasings discussed in the pre- 
vious section in terms of GW detection. In order to do 
so, we calculate overlaps and fitting factors [28l . |43 | be- 
tween point-particle waveforms, which we treat as our 
templates, and "real" signals, i.e. waveforms in which 
we include tidal effects to the best of our knowledge by 
means of Eq. ^T\\ . Our templates and signals thus have 
the form 



Ae 



(62) 
(63) 



where A denotes the GW amplitude and V'tj '0t^^, and 
SipT are the tidal contributions to the Fourier phase given 
in Eqns. ([SS)) . ([57]) . and (pO]) . respectively. 

Given two signals hi and ft,2, their noise- weighted inner 
product is 



where S'n(/) is the power spectral density of a given 
detector, the lower integration bound /start depends on 
the detector one is examining, and for the upper bound 
/cut we use twice the orbital frequency /™^ discussed 
in the previous section. More specifically, we consider 
second and third generation detectors, such as Advanced 
Virgo/LIGO, and the Einstein Telescope, setting /start to 
20 Hz and 10 Hz, respectively, and follow for S'n(/). 
The inner product allows us to define the overlap between 
two signals, that is, their normalized inner product, max- 
imized over time and phase shifts: 



0[hi,h2 



max 



{hi\h2 



(65) 



(64) 



where t and (p are the time and phase offsets between 
the two waveforms. This quantity is useful, for example, 
when describing in quantitative terms the effectiveness 
of a waveform model in detecting a physical waveform. 
Furthermore by maximising 0[hi,h2] over the intrinsic, 
physical parameters of, say, hi, one obtains the fitting 
factor. In our case, the point-particle waveform tem- 
plates, hpp, depend on the binary total mass and sym- 
metric mass ratio, so that we denote the fitting factor 
with !F!F{m, v). 

As discussed in [2^, when calculating the number of 
missed events in a search performed with a discrete tem- 
plate bank, one must take into account the template 
bank spacing and the fitting factor TT(m, v) between 
the template GW model and the real signal model. In 
the case of LIGO/ Virgo template banks, this may be 
done by subtracting a maximum mismatch of 0.03 to the 
fitting factors we calculate: this yields the effective fit- 
ting factors ETJF(m, v) = TT(m, v) — 0.03 which must 
then be used to determine the fraction of missed events 
1 — £TT{m,v)^ . [23 additionally discusses a criterion 
for template waveform accuracy: 1 — TT(m, v) must be 
smaller than 0.005, so that 1 — f J-"J-"(m, v) < 0.035 for 
LIGO/ Virgo. 

Our results may be summarized as follows: 

• We find that the lowest overlaps occur in the case 
of low mass ratios as already noted in ; this hap- 
pens for two reasons, in general: because low total 
masses and mass ratios enhance the dephasing in 
Eqs. ([57]) and (pH)) and because /cut is higher 
for lower mass systems, thus allowing the phase dif- 
ference originating from tidal distortions to accu- 
mulate over a broader frequency range. Further- 
more, in the q = 1 case the phase difference orig- 
inates from the distortion of two NSs instead of a 
single one. 

• For BH-NS binaries, the 0[hpp,hT,s] overlaps are 
> 0.997 for second generation detectors and > 
0.995 for the Einstein Telescope, as already found 
in 0. This means that (1) point-particle tem- 
plates are within the needed waveform accuracy, 
since l-J"J"<l-e'< 0.005, and that (2) the 
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total number of missed events including the effect 
of template bank spacing, is 1 — £FF{rn^ vY < 
1 - (O - 0.03)3 < 10%. 

For second generation detectors and NS-NS bi- 
naries, we find that, depending on the EOS, 
the overlap C[ft.pp, ft.T,(5] varies within the range 
[0.963,0.997] for ms ' = 1.2Mq and within 
[0.982, 0.999] for = IAMq; the lower hmits cor- 
respond to larger deformabilities, i.e. to the stiffer, 
PS EOS. For third generation detectors, the inter- 
vals are, instead, [0.943,0.992] and [0.969,0.997]. 

By computing the fitting factors for NS binaries 
with the PS EOS, we obtain FF{m,v) = 0.991 
for 7712 = 1.2Af0 and FF{m,v) = 0.996 for m2 = 
1.4M0, in the case of second generation detec- 
tors. For third generation detectors, one has, in- 
stead, FF{rn,v) = 0.985 for m2 = 1.2A/0 and 
J"J"(m, v) = 0.992 for ma = 1.4Afo. 

In the case of Advanced LIGO /Virgo and of a par- 
ticularly stiff EOS for matter in the NS interior, the 
total number of missed events corresponding to the 
above fitting factors, 1 — £J-T{m, v)^, may thus be 
as high as 11% for low mass NS binaries, whereas it 
would roughly be 10% for a canonical m = 2.8Mq 
equal mass NS binary. The Einstein Telescope, on 
the other hand, would miss up to ~ 13% low mass 
NS-NS inspirals and up to 11% canonical NS inspi- 
rals. Our results are summarized in Tables HIl HID 

If the EOS of NS matter is very stiff, point-particle 
templates of binary NS inspirals would not meet 
the required accuracy for Advanced Virgo/LIGO 
in the case of a low total mass, since the fitting 
factor we obtain yields 1 — J"J"(to, v) > 0.005 [28j . 
The same is true for equal mass NS inspirals and 
the Einstein Telescope, on the high stiffness end of 
possible EOSs. As previously mentioned, the PS 
EOS is an extreme case of stiff EOS. Nevertheless, 
one should consider this worst (in terms of missed 
events and waveform template accuracy) case sce- 
nario and establish if, and eventually how, we may 
remedy; this is important also because a lot of in- 
teresting physics is associated with low mass NS 
binary mergers, e.g. [i^ H^- If one imagines to 
use templates /it.pn that include tidal corrections 
up to (relative) IPN order in the Fourier phase by 
means of the Love Number, i.e. 

^T.iPN = y^e^'^'PP+^'^T+^T'"' (66) 

all overlaps ©[/it.ipn, ^t,i5] would differ from unity 
by less than a part in one-thousand, both for second 
and third generation detectors. To the best of our 
knowledge in the modelling of the real signal, the 
inclusion of Love Number dependent tidal terms in 
GW templates would therefore greatly help in con- 
straining the NS EOS in post-processing analysis, 



confirming the result of . Building gravitational 
waveforms within the adiabatic approximation (i.e., 
including if;^ and tp^^) is thus very reliable. 

• Let us now suppose that the EOS of NS matter is 
stiff and that an to = 2AMq equal mass binary NS 
inspirals close enough to Earth. The detection of 
its gravitational radiation by a second generation 
detector with an imaginarily infinitesimally spaced 
template bank would hit a point-particle template 
with total mass to = 2.446Mq and symmetric mass 
ratio V = 0.242, since these are the values that 
maximise the overlap ©[/ipp, h^^s]'- this means that 
tidal effects contribute with a 2% error to the to- 
tal mass measurement and with a 3% error on the 
symmetric mass ratio measurement. In the case 
of an TO = 2.8Mq equal mass inspiral, the tidal 
contributions to the errors on the total mass and 
the symmetric mass ratio measurements would be 
1% and 2%, respectively. The same values hold for 
third generation detectors. 



EOS 


m(M0) 


0[hpp,hT,5\ 


TT(m, v) 


missed 


APR 


2.4 


0.997 


1.000 


0% 




2.8 


0.999 


1.000 


0% 


PS 


2.4 


0.963 


0.991 


11% 




2.8 


0.982 


0.996 


9.8% 



TABLE II. We report the values of overlaps 0\}ipp, hT,s] and 
fitting factors J-J-{m, u) between point-particle templates and 
waveforms in which we include tidal effects evaluated as in 
Eq. (|6ip . We also show the total number of missing events. 
All values refer to NS-NS systems and second generation de- 
tectors AdVirgo/LIGO. 



EOS 


m(M0) 


0[hpp,hT.6] 


TJ'{m, u) 


^ missed 


APR 


2.4 


0.992 


1.000 


0% 




2.8 


0.997 


1.000 


0% 


PS 


2.4 


0.943 


0.985 


13% 




2.8 


0.969 


0.992 


11% 



TABLE III. The same quantites of Table |Ill but for third 
generation detector ET. 



VI. CONCLUDING REMARKS 

In this paper we improve the PNA model recently de- 
veloped in [l6| , by computing the tidal tensor of a binary 
system up to 2 PN order. Using this approach, we study 
the dynamical evolution of NS deformations in compact 
binary inspirals, finding a number of interesting results. 

We find that the ratio between the quadrupole tensor 
qij and the tidal tensor increases during the last stages 
of the inspiral, and we show that the Love number has to 
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be considered as the asymptotic value of this ratio for suf- 
ficiently large orbital separation. This increase, however, 
only marginally affects the phase of the emitted signal, 
thus assessing the validity of the adiabatic approximation 
when modeling the gravitational wave signal. 

We provide a semi-analytical expression of the Love 
number, in terms of the scalar quadrupole moment and 
of the integral of pressure over the stellar volume, both 
computed for the star at isolation, i.e. for a spherically 
symmetric configuration (Eq. (I49|) ). This expression pro- 
vides a physical insight on the Love number and a for- 
mulation alternative to that given in [l^ . only in terms 
of quantities which refer to the structure of the star in 
its spherical equilibrium configuration. 

In addition, we estimate the reliability of point par- 



ticle templates with respect to a "real" signal modeled 
by means of our approach, finding that such templates 
marginally fail to meet the standard accuracy require- 
ments (1 - FF{m, v) < 0.005) for NS-NS inspirals with 
very stiff equation of state. 
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